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" Sufvi de la deformation d'une structure lineique definie sur une image d'une 
sequence d'images d'un organe deformable dans le temps. " 

Description : 

5 

invention concerne un procede de traitement d'images appartenant a une sequence d'au 
moins deux images presentant une surface representative d'un organe ou d'une partie 
d'organe deformable dans le temps et appelee surface d'organe, ladite surface incluant des 
points caracteristiques, notes points de marquage, qui se correspondent d'une image a une 
10 autre de la sequence. L' invention concerne aussi un appareillage de traitement d'images 
pour mettre en oeuvre le procede decrit ci-dessus. 

L'invention trouve son application dans le domaine du traitement d'images medicales. Le 
procede est en particulier applicable sur les images d'organe marque par modulation spatiale 
de resonance magnetique. Ce marquage est visible sur les images sous la forme de lignes de 
15 marquage ayant des points d'intersection. Les lignes de marquage se deferment en suivant 
la deformation de I'organe. Lesdits points d'intersection sont alors choisis comme points de 
marquage car les correspondences de ces points d'une image a I'autre sont aisement 
detectables. 

20 Un procede de traitement d'images d'organe marque par modulation magnetique en vue de 
quantifier la deformation de I'organe est deja connu de Petat de la technique par la 
publication de Matthias Stuber et al. intitulee "Quantification of the local heartwall motion^by 
magnetic resonance myocardial tagging". Dans ce document, les points d'intersection entre 
les lignes de marquage sont determines par des algorithmes de "snakes". Par ailleurs, ce * 

25 procede met en ceuvre une moyenne des angles par rapport au centre de gravite pour 
quantifier la rotation et la contraction. 

Ce procede presente des inconvenients. D'abord, 1'utilisation des algorithmes de snakes pour 
la determination des points d'intersection entre les lignes de marquage manque de precision. 

30 Ensuite, la mise en ceuvre de la moyenne des angles par rapport au centre de gravite pour 
quantifier la rotation et la contraction resulte aussi en un resultat imprecis et 
mathematiquement incorrect. Par consequent, des extrapolations ulterieures ne sont pas 
fiables pour les points autres que les points de marquage. Or, ces extrapolations sont 
necessaires lorsque I'on souhaite suivre le mouvement d'une structure lineique particuliere 

35 qui est definie sur une Image et qui est amenee a se deformer avec I'organe dans la 
sequence d'images. 
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* Suivi de la deformation d'une structure lineique definie sur une image d'une 
sequence d'images d'un organe deformable dans le temps. " 

Description : 

Llnvention concerne un procede de traitement d'images appartenant a une sequence d'au 
moins deux images presentant une surface representative d'un organe ou d'une partie 
d'organe deformable dans le temps et appelee surface d'organe, ladite surface incluant des 
points caracteristiques, notes points de marquage, qui se correspondent d'une image a une 
autre de la sequence. L 1 invention concerne aussi un appareillage de traitement d'images 
pour mettre en oeuvre le procede decrit ci-dessus. 

invention trouve son application dans le domaine du traitement damages medicales. Le 
procede est en particulier applicable sur les images d'organe marque par modulation spatiale 
de resonance magnetique. Ce marquage est visible sur les images sous la forme de lignes de 
marquage ayant des points d'intersection. Les lignes de marquage se determent en suivant 
la deformation de I'organe. Lesdits points d'intersection sont alors choisis comme points de 
marquage car les correspondences de ces points d'une image a I'autre sont aisement 
detectables. 

Un procede de traitement d'images d'organe marque par modulation magnetique en vue de 
quantifier la deformation de I'organe est deja connu de I'etat de la technique par la 
publication de Matthias Stuber et al. intitulee "Quantification of the local heartwall motion by 
magnetic resonance myocardial tagging". Dans ce document, les points d'intersection entre 
les lignes de marquage sont determines par des algorithmes d'approximation bases sur un 
calcul de potentiel (en anglais, « snakes »). Par ailleurs, ce procede met en oeuvre une 
moyenne des angles par rapport au centre de gravite pour quantifier la rotation et la 
contraction. 

Ce procede presente des inconvenients. D'abord, Tutilisation des algorithmes 
d'approximation bases sur un calcul de potentiel pour la determination des points 
d'intersection entre les lignes de marquage manque de precision. Ensuite, la mise en oeuvre 
de la moyenne des angles par rapport au centre de gravite pour quantifier la rotation et la 
contraction resulte aussi en un resultat imprecis et mathematiquement incorrect. Par 
consequent, des extrapolations ulterieures ne sont pas fiables pour les points autres que les 
points de marquage. Or, ces extrapolations sont necessaires lorsque Ton souhaite suivre le 
mouvement d'une structure lineique particuliere qui est definie sur une image et qui est 
amenee a se deformer avec I'organe dans la sequence d'images. 
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Un but de ('invention est de fournir un precede pour etudier la deformation d'un organe a 
partir des donnees de (Image sur une structure lineique definie sur une image de la 
sequence. Un but particulier de 1'invention est de fournir un procede pour suivre, avec 
fiabilite, la deformation d'une telle structure lineique definie sur une image. Le suivi d'une 
structure lineique est en particulier interessant pour etudier le mouvement des parois de 
I'organe qui apparaissent sur llmage sous la forme d'un contour. Llnvention s'applique 
particulierement au cceur qui presente deux parois majeures visibles sur les images, par 
exemple obtenues par resonance magnetique : I'endocarde et I'epicarde. 
Selon 1'invention, les buts cites sont atteints, et les problemes presentes par I'etat de la 
technique sont resolus, par un procede conforme au paragraphe introductif et caracterise en 
ce qull comprend des etapes de: 

- definition d'une structure lineique sur une des images de la sequence, 

- calcul des positions des points de marquage sur au moins deux images, successives ou 
non, 

- determination des parametres d'une expression mathematique explicite de la deformation 
de i'organe ou de ia partie de I'organe observee entre les deux images, a partir des positions 
d'un ensemble de points de marquage sur les deux images et des positions des points de la 
structure lineique, 

- application de ladite expression mathematique explicite a la structure lineique pour definir 
la forme de la structure lineique apres deformation de I'organe entre les deux images. 

Le procede selon Pinvention presente de nombreux avantages parmi lesquels les avantages 
de la regularisation par une expression mathematique explicite pour quantifier la 
deformation. Cette regularisation remplace I'approche par interpolation du mouvement en ' 
chaque point qui est le procede choisi dans le document de I'etat de la technique cite. x 
L'expression mathematique obtenue est alors au moins valable en tout point de ia structure 
lineique definie lors de la premiere etape. Cette approche estime un mouvement au sens 
tres rigoureux des moindres carres avec une regularisation mathematique explicite. Cette 
regularisation corrige le bruit. 

Dans une mise en oeuvre avantageuse de llnvention, la structure lineique realise une 
segmentation de Umage, ladite segmentation etant ainsi suivie d'une image a Tautre. La 
structure lineique definit plusieurs zones dans Umage et les calculs de la. deformation de 
I'organe peuvent alors etre realises uniquement sur une zone dlnteret particulier. Par 
exemple, dans le cas du cceur, la definition de deux structures lineiques suivant 
respectivement I'endocarde et Tepicarde permet de ne realiser les calculs que sur la surface 
de Umage reeilement representative de I'organe. Cela diminue les temps de calcul et fiabilise 
le resultat car les mouvements "parasites" des structures n'appartenant pas a I'organe ne 
sont pas pris en compte. 
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Dans un premier mode de realisation de finvention, I'expression mathematique de la • 
deformation de I'organe ou de la partie de I'organe observee entre les deux images est 
determinee pour une surface incluant la structure lineique a partir des positions sur les deux 
images d'un ensemble de points de marquage qui contient au moins les points de marquage 
5 presents sur ladite surface. La structure lineique peut etre incluse strictement dans la 
surface ou etre le contour de cette surface. 

Dans un second mode de realisation, I'expression mathematique de la deformation de 
I'organe ou de la partie de I'organe observee entre les deux images est determinee 
uniquement pour les points de la structure lineique a partir des positions sur les deux images 
10 d'un ensemble de points de marquage qui sont ponderes en fonction de leur eloignement 

par rapport a la structure lineique. Les poids attribues aux points de marquage sont d'autant 
plus importants que le point est proche de la structure lineique. Ces poids, utilises dans une 
approximation aux moindres carres, permettent d'obtenir une expression mathematique tres 
fiable pour les points de la structure lineique ou elle est definie. 

15 

Llnvention sera mieux comprise a la lumiere de la description suivante de quelques modes 
de realisation, faite a titre d'exemples et en regard des dessins annexes, dans lesquels :■■ 

la figure 1 represente un diagramme d'un procede de traitement dlmages selon 

Invention, 

20 - la figure 2 presente une image CSPAMM d'un coeur, ladite image appartenant a une 

sequence dlmages et etant prise apres llmpulsion de magnetisation, 

la figure 3 illustre un premier mode de realisation de finvention, 

la figure 4 illustre un second mode de realisation de Invention, 

la figure 5 represente un appareillage de traitement dlmages selon une mise en oeuvre 
25 particuliere de Hnvention, 

- les figures 5a, 6b, 6c et 6d represented la deformation d'une structure lineique 

respectivement sans correction du mouvement global, avec correction de translation 

globale, avec correction de translation et de rotation globales, avec correction de 

translation, de rotation et de contraction globales, 
30 - la figure 7 represente un appareil de saisie et de traitement dlmages selon nnvention. 



La figure 1 represente un diagramme d'un procede de traitement dlmages selon (Invention. 
Ledit procede est applicable sur des images appartenant a une sequence dlmages d'au 
moins deux images IM(tl) et IM(t2) prises a deux instants tl et t2 d'un organe ou d'une 
35 partie d'un organe amene a se deformer dans le temps- Ledit organe ou ladite partie de 
I'organe est visible sur les images sous la forme d'une surface appelee surface de I'organe. 
La surface de I'organe inclut des points caracteristiques dont les correspondences sont 



1 er depot 




4 

determinees d'une image a une autre de la sequence. Lesdits points caracteristiques sont 
notes points de marquage. L'organe ou la partie d'organe peut, par exemple, etre marque(e) 
par modulation spatiale de resonance magnetique. Dans la suite, llnvention est decrite plus 
particulierement dans le cas de ce marquage par modulation spatiale de resonance 
magnetique. 

La technique de marquage par modulation spatiale de resonance magnetique inclut en 
particulier les techniques SPAMM et CSPAMM permettant d'obtenir des images ou le 
marquage est visible sur les images suivant des lignes de marquage qui peuvent etre de 
differentes geometries lorsqu'dles sont generees dans I'organe (droites, courbes...). Lesdites 
lignes de marquage se determent en suivant la deformation materielle de I'organe. Sur les 
images obtenues grace aux techniques evoquees ci-dessus, les lignes correspondant aux 
minima de magnetisation spatiale sont des lignes sombres et sont aisement reperables. 
La modulation spatiale de resonance magnetique est en general utilisee en prenant des 
series dlmages de I'organe marque a des temps successifs et reguliers. Ces series dlmages 
sont denommees sequences dlmages et la deformation de I'organe est observee grace a la 
deformation des lignes de marquage qui constituent une sorte de trame attachee a I'organe. 
Ladite trame peut avoir divers aspects : droites paralleles, grille constitute de droites suivant 
deux directions... Une technique denommee "Slice-Following" permet de suivre la 
deformation d'une section de I'organe meme si le plan de cette section se deplace dans une 
direction substantiellement perpendiculaire a ce plan pendant la sequence. 
Dans le cas d'une deformation periodique de I'organe, plusieurs sequences dlmages du . 
meme organe prises pour des deformations similaires successives montrent la meme 
deformation sur chaque image similaire, e'est-a-dire, echantillonnee au meme instant au sein 
de la periode de deformation. Dans ce cas, les images similaires peuvent etre combinees de. 
maniere a ce que les trames des deux images soient visibles sur la nouvelle image resultant 
de la combinaison. On definit ainsi une nouvelle sequence dlmages en effectuant cette 
combinaison sur I'ensemble des images de la sequence. Cette nouvelle sequence contient en 
general plus dfnformation que les sequences d'origine plus simples. 
En fonction des temps d'acquisition des sequences qui dependent du marquage choisi, il 
peut etre avantageux, soit de travailler sur une seule sequence dlmages de I'organe marque 
avec un marquage complexe, soit de travailler sur une combinaison de plusieurs sequences 
dlmages (generalement deux) de I'organe chaque fois marque avec un marquage simple, 
ladite combinaison definissant une nouvelle sequence utilisee ensuite dans le procede de 
traitement dlmages selon llnvention. 

La figure 2 presente une image d'un coeur, ladite image appartenant a une sequence 
dlmages et etant prise environ 9 ms apres llmpulsion de magnetisation. Deux ensembles de 
lignes paralleles correspondant a des minima dlntensite lumineuse sont observes, les lignes 
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paralleles d'un ensemble etant perpendiculaires aux lignes paralleles de I'autre ensemble. La 
sequence dlmages dont est tire cet exemple presente ainsi des lignes de marquage suivant 
deux directions distinctes et peut, ainsi, etre le resultat, sort d'une saisie directe d'une 
sequence damages de I'organe marque suivant les deux directions, soit de la combinaison de 
5 deux saisies de sequences, chacune des deux sequences etant marquee suivant une des 
deux directions. 

Deux types de lignes de marquage peuvent etre utilisees dans un procede suivant 
llnvention. Sur la figure 2, les lignes de marquage correspondant a des minima dlntensite, 
c'est-a-dire a des minima de resonance magnetique sont facilement reperables. Les lignes 

10 de marquage correspondant a des maxima dlntensite et correspondant a des maxima de 
resonance magnetique sont cependant egalement detectables, meme si leur detection est 
moins aisee. Par exemple, par derivation du profil dlntensite de llmage i( est possible de 
localiser les lignes correspondant aux maxima de magnetisation. L'utilisation de ces deux 
types de lignes de marquage augmente Information sur I'image car la trame du marquage 

15 est plus serree : le nombre de points dlntersection entre lignes de marquage et done le 
nombre de points de marquage est plus important 

A propos de la qualite des prises dlmages d'une sequence, la technique CSPAMM permet 
notamment d'obtenir un contraste persistant sur une sequence. Cela est en particulier utile 
dans le cas ou les lignes de marquage correspondant aux maxima de resonance magnetique 

20 sont utilisees, le contraste persistant aidant au reperage des minimums dlntensite. 

Le procede de traitement dlmages suivant llnvention traite les images ou sont presents des 
points de marquage pour lesquete il est possible d'etablir des correspondances d'une image 
a I'autre. Dans le cas de la modulation spatiale de resonance magnetique, les lignes de 
marquage sont telles qull existe des points dlntersection entre plusieurs lignes de 

25 marquage. II est aise d'etablir des correspondances d'une image a I'autre pour ces points 
d'intersection qui sont nommes par la suite points de marquage MP. lis peuvent etre des 
points d'intersection entre lignes de marquage de toute forme et peuvent etre directement 
visibles sur la trame de marquage ou n'etre visible qu'apres une combinaison de plusieurs 
sequences, donnant une nouvelle sequence dlmages. 

30 Le procede seion llnvention inclut une etape DEF de definition sur une image IM(ti) d'une 
structure lineique LS(ti). Dans I'application du procede dans un appareillage de traitement 
d'images suivant llnvention, cette etape peut par exemple etre realisee manuellement par 
I'utilisateur ou automatiquement sur une image quelconque de la sequence. Cette image 
IM(tl) peut etre choisie a nlmporte quel instant tl de la sequence car les deformations 

35 observees d'une image a Tautre peuvent etre, une fois determinees, appliquees 

indifferemment dans le sens chronologique ou dans le sens contraire par simple inversion 
des deformations. Dans le cas ou la structure lineique permet de suivre la deformation d'une 
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paroi de I'organe observe, I'image sera avantageusement choisie comme etant celle ou ladite 
paroi est la plus visible : le praticien peut, dans ce cas, definir la structure lineique en 
suivant simplement visuellement le contour sur I'image. 

Le procede selon llnvention inclut une etape de calcul CALC des positions des points de 
marquage. Cela peut, dans I'exemple de la modulation spatiale de resonance magnetique, 
etre realise par utilisation du procede decrit dans la demande de brevet incluse par 
reference : No 00403028.4 du 31 octobre 2000. Dans cette demande de brevet des points 
candidate pour appartenir a une ligne de marquage donnee sont detectes avant que des 
moyens de prevision du mouvement de la ligne de marquage soient utilises pour identifier la 
ligne de marquage et les points lui appartenant et qu'une equation de la ligne ne soit 
calculee. Une fois les equations des lignes de marquage determinees, le calcul des positions 
des points de marquage est aise. Les points de marquage MP(tl) et MP(t2) sont les points 
pour lesquels une correspondance est etablie d'une image IM(tl) a llmage IM(t2), les deux 
instants tl et t2 pouvant etre ou non successife dans la sequence dlmages et tl pouvant 
etre posterieur ou anterieur a t2. 

Le procede selon llnvention inclut ensuite une etape de determination DET d'une expression 
mathematique f de la deformation de I'organe ou de la partie de I'organe observee entre 
llmage IM(tl) et IM(t2) a partir d'un ensemble MP' de points de marquage. 
Un premier mode de realisation de I'etape de determination DET est illustree sur la figure 3. 
La figure 3 represente une image IM sur laquelle les positions des points de marquage MP 
sont calcules. L'expression mathematique f de la deformation de I'organe ou de la partie de: 
I'organe observee entre cette image et une seconde image est determinee pour tous points \ 
d'une surface S incluant la structure lineique LS a partir des positions sur les deux images ?? 
d'un ensemble de points de marquage MP 1 , (edit ensemble de points de marquage MP' etant 
dans cet exemple, Tensemble des points de marquage presents sur ladite surface S. La 
structure LS etant incluse dans la surface S, la deformation de la structure lineique LS de 
cette image a la seconde est obtenue en appliquant l'expression mathematique f a LS. 
Dans ce premier mode de realisation, l'expression mathematique est determinee sur toute la 
surfaces. 

Un deuxieme mode de realisation, plus precis, de I'etape de determination DET est illustree 
sur la figure 4. La figure 4 represente une image IM sur laquelle les postions des points de 
marquage MP sont calcules. L'expression mathematique de la deformation f de I'organe ou 
de fa partie de I'organe observee entre les deux images est determinee uniquement pour les 
points de la structure lineique LS a partir des positions sur les deux images d'un ensemble 
de points de marquage MP 1 , lesdits points de marquage MP 1 etant une. partie des points de 
marquage MP ou la totalite de ces points de marquage MP ponderes en fonction de leur 
eloignement par rapport a la structure lineique LS. Le poids a, attribue a un point i 
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appartenant a I'ensemble des points de marquage MP* est d'autant plus important que le 
point de marquage est proche de la structure lineique LS. Par exemple le poids a, peut etre 
de la forme : 

1 

a i = 1 vT - 

[distancef point i;LS )) 

Les parametres de I'expression mathematique f sont generalement determines par 
approximation aux moindres carres, en utilisant ou non une ponderation des poids selon le 
mode de realisation. Dans le cas de points ponderes, I'approximation aux moindres carres 
revient a chercher la fonction f realisant : 

f ^ 
Inf ^ai(f(pointi(tl))~pointi(t2)) 2 
\MP f 

Cette approximation est realisee a partir des positions des points de marquage dont les 
positions sont connues sur les deux images et done, dont le deplacement est connu entre tl 
et t2, en se referant a nouveau a la figure 1. En utilisant une ponderation des points de 
marquage MP 1 , la regularisation sera ainsi meilleure au voisinage de la structure lineique. 
[/expression mathematique peut, par exemple, etre une similitude qui rend compte diss 
deformations rigides. 

Dans un mode de realisation particulier, Texpression mathematique de la deformation est 
definie dans le plan complexe. La deformation peut etre definie dans un repere polaire. Dans 
le cas particulier ou I'organe est un coeur, il est aise de definir un centre positionne 
approximativement au centre de gravite de la surface de I'organe telle que vue sur llmage. 
Dans le mode de realisation prefere de llnvention, I'expression mathematique est derivee 
d r une similitude que Ton peut ecrire sous une forme relative a un point achoisi de maniere 
quelconque, mais generalement choisi comme etant approximativement le bar/centre de la 
surface observee : 

f(z) = | z-o I (W V d , 0 = arg(z-o). 
L'expression de la similitude est modifiee en introduisant une serie de Fourier dans 
I'expression pour rendre compte de deformations plus elastiques qu'une simple similitude : 

N 

+ d , 0=arg(z-o), (a k )eC 2N ^ 



f(z) = \z-o\ 



a k e 



k=-N 
K k*Q 



Cette expression rend compte d'une deformation semi-elastique globale dans un repere 
centre. Elle est definie par 2x(2N+l)+2 parametres reels qui sont definis a partir des points 
de marquage de I'ensemble MP'. Ces parametres sont en general surevalues car il y a plus 
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de points de marquage dans MP* que de parametres. C'est en particulier le cas lorsque les 
lignes de marquage correspondant aux maxima de resonance magnetique sont utilisees, le 
nombre de points dlntersection entre les lignes etant alors important. La surevaluation des 
parametres permet de lisser le bruit. 

Dans le cas du coeur, 1'endocarde est notablement plus contractile que I'epicarde et par 
consequent, plus on approche du centre du myocarde, plus Tamplitude du mouvement radial 
est importante. Un terme correctif fonction du rayon et Tangle polatre est avantageusement 
ajoute a posteriori a la determination de 1'expression mathematique explidte f pour ajouter 
une dependance en rayon. Ce terme correctif est lui aussi determine en utilisant les 
positions des points de marquage MP' sur les deux images mais cette determination est 
realisee apres la determination des parametres a k de la deformation f(z). 
Dans une mise en oeuvre avantageuse, le terme correctif est defini par secteurs angulaires s 
de Hmage et est de la forme : 



Y (r,0)^ 



r,(rh 



Le terme y s (r) est un polynome en r independant de 9 defini sur le secteur angulaire s de, 
I'image en fonction des positions des points de marquage sur le secteur angulaire s sur les 
deux images. Llnterpolateur polynomial de Lagrange est ensuite utilise pour tenir compte de 
la. dependance en G, 0 S etant I'angle central du secteur angulaire s. Si le terme correctif etait 
determine en meme temps que la fonction f il existerait une multiplicity d'ecriture de 
I'ensemble sauf si la fonction f etait contrainte. Or f donne une expression mathematique de 
la deformation globale. Le terme correctif represente le comportement physiologique et non 
pas la deformation cinematique. Globalement, I'elasticite de I'expression mathematique est 
controlee par le choix de N, ordre de Fourier (N=3 est generalement suffisant) et, dans le 
cas ou un terme correctif est introduit, par le nombre de secteurs angulaires consideres pour 
I'approche du terme correctif et le choix du degre du polynome y s . 

Dans le second mode de realisation de Tetape de determination, il est possible d'avoir une 
determination tres precise de la deformation en appliquant dans un premier temps une 
etape de determination d'une expression mathematique suivant le premier mode de 
realisation, c'est-a-dire sur une surface, en utilisant un terme correctif. Ce terme correctif, 
une fois determine, est ensuite utilise dans I'etape de determination de I'expression 
mathematique determinee uniquement sur la structure lineique. Les deux modes de 
realisation peuvent ainsi etre mis en oeuvre independamment ou etre combines pour obtenir 
des donnees complementaires ou une meilleure precision. 

La figure 5 represente un appareillage DEV de traitement d'images selon une mise en oeuvre 
particuliere de nnvention. Cet appareillage DEV est en relation avec des moyens 
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d'acquisition ACQ de sequences SIM de X images. Cet appareillage DEV inclut des moyens 
DEF de definition d'une structure lineique LS(t|) sur une image IM(tj). Cet appareillage DEV 
inclut des moyens de calcul CALC des positions des points de marquage sur deux images 
successives ou non, et des moyens de determination DET des parametres d'une expression 
mathematique explicite de la deformation de I'organe ou de la partie de I'organe observee 
entre les deux images a partir des positions d'un ensemble de points de marquage sur les 
deux images, ladite expression mathematique etant au moins determinee pour les points de 
la structure lineique. Dans la mise en oeuvre particuliere representee sur la figure 5, le 
procede decrit dans la figure 1 pour deux images successives ou non est itere sur I'ensemble 
des images successives IMft) de la sequence damage SIM. 

Comme la structure lineique est definie sur une image IM(tj) de la sequence qui est le plus 
souvent differente de IM(to), le traitement des images ne peut etre fait a partir de tb 
directement. II est necessaire de reordonner la sequence SIM de maniere a commencer a 
traiter llmage ou est definie la structure lineique en premier, quitte a ensuite retablir I'ordre 
d'origine pour la presentation des resultats. Les deformations peuvent en effet etre calculees 
entre deux images quelconques. L'utilisation d'une memoire MEM permet aussi de stocker 
les expressions des deformations pour ensuite les reutiliser. Par exemple, il est possible de 
combiner entre elles les expressions des deformations dlmages successives pour trouver la 
deformation entre deux images non successives. 

Apres une initialisation pour un compteur initialise a i=0, (Image IM(tj=to) est extraite de la 
sequence SIM. Les positions des points de marquage MP(tj=to) sont calculees et stockees 
dans la memoire MEM, le processus decrit ci-apres est amorce avec i=l. Cette initialisation 
n'est pas explicitement representee sur la figure car elle est le cas particulier du schema 
general ou 1=0, ou f est fonction nulle et MP(t hl )=MP(t j =tb). Une image IMft) est extraite de 
la sequence dlmages SIM a priori reordonnee. Les positions des points de marquage MP(tj) 
sont calculees par des moyens de calcul CALC. Ces positions sont stockees dans une 
memoire MEM et sont fournies a des moyens de determination DET d'une expression 
mathematique de la deformation. Les positions d'un ensemble de points de marquage 
MP'(thi) de llmage precedente IM(tn) sont extraits de la memoire MEM et fournis aux 
moyens DET pour la determination de I'expression de la deformation f(Wtj). 
Dans la mise en oeuvre representee sur la figure 5, Texpression mathematique de la 
deformation est ensuite stockee dans la memoire MEM et le compteur est increments a 
i=i+l. La deformation f(Wtj) est utilisee avec LSfti-O dans des moyens duplication AP de 
I'expression mathematique de fa deformation a ISO+d pour determiner LSft) qui est stockee 
dans la memoire MEM et peut etre representee sur llmage IMft). 
La figure 6 illustre un dernier perfection nement du suivi de la deformation d'une structure 
lineique definie sur une image incluant des points de marquage. Dans une mise en oeuvre 
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avantageuse, les mouvements globaux de la structure lineique sont en effet calcules puis 
extraits de I'expression appliquee a la structure lineique. La figure 6 represents la 
deformation de la structure lineique en representant deux structures lineiques : une 
premiere structure LS(tl) correspondant a la structure sur une image donnee et une 
deuxieme structure LS(t2) qui a subi la deformation de I'organe entre les instants tl et t2. Le 
deplacement specifique de douze points est represente par des lignes, telle la ligne L reliant 
les points correspondents d'une structure a I'autre. Sur la figure 6a, aucune correction n'est 
effectuee, sur la figure 6b, est appliquee une correction de la translation globale de la 
structure lineique, sur la figure 6c, est appliquee en plus de la correction de translation, une 
correction de la rotation globale, sur la figure 6d, est appliquee en plus une correction de 
contraction. La figure resultante donne des informations sur la deformation locale 
particuliere de I'organe, ce sont ces informations qui permettent de detecter d'eventuelles 
anomalies dans le comportement de I'organe. 

La figure 7 represente un appareil APP de saisie dlmages, ledit appareil comprenant des 
moyens ACQ d'acquisition de sequences d'au moins deux images d'un organe ou d'une 
partie d'un organe amene a se deformer dans le temps, ledit organe ou ladite partie de 
I'organe etant visible sur les images sous la forme d'une surface appelee surface de I'organe, 
ladite surface incluant des points caracteristiques dont les correspondances sont 
determinees d'une image a une autre de la sequence, des moyens REP de representation 
visuelle de ces images pouvant comprendre un mode video pour suivre la deformation au 
cours de la sequence, un appareillage DEV de traitement dlmages tel que decrit auparavant 
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Revendications : 

1. Procede de traitement d'images appartenant a une sequence d'au moins deux images 
presentant une surface representative d'un organe ou d'une partie d'organe deformable 
dans le temps et appelee surface d'organe, ladite surface incluant des points 
caracteristiques, notes points de marquage, qui se correspondent d'une image a une autre 
de la sequence, ledit procede comprenant des etapes de : 

~ definition d'une structure lineique sur une des images de la sequence, 
_ calcul des positions des points de marquage sur au moins deux images, successives 
ou non, 

- determination des parametres d'une expression mathematique explicite de la 
deformation de I'organe ou de la partie de Torgane observee entre les deux images, a partir 
des positions d'un ensemble de points de marquage sur les deux images et des positions des 
points de la structure lineique, 

- application de ladite expression mathematique explicite a la structure lineique pour 
definir la forme de la structure lineique apres deformation de I'organe entre les deux images. 

2. Procede de traitement d'images selon la revendication 1, caracterise en ce que ledit 
organe est marque par modulation spatiale de resonance magnetique, ledit marquage etant 
visible sur les images sous la forme de lignes de marquage, lesdites lignes de marquage se 
deformant en suivant la deformation de I'organe et etant telles qull existe des points 
dlntersection entre lesdites lignes de marquage, lesdits points dlntersection etant les points 
de marquage. 

3. Procede de traitement d'images selon I'une des revendications 1 ou 2, caracterise en 
ce que la structure lineique definit une segmentation de 1'image, ladite segmentation etant 
ainsi suivie d'une image a I'autre. 

4. Procede de traitement d'images selon I'une des revendications 1 a 3, caracterise en ce 
que ('expression mathematique de la deformation de I'organe ou de la partie de I'organe 
observee entre les deux images est determinee pour une surface incluant la structure 
lineique a partir des positions sur les deux images d'un ensemble de points de marquage, 
ledit ensemble de points de marquage contenant au moins les points de marquage presents 
sur ladite surface. 

5. Procede de traitement d'images selon I'une des revendications 1 a 3, caracterise en ce 
que I'expression mathematique de la deformation de I'organe ou de la partie de I'organe 
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observee entre les deux images est determinee uniquement poiir les points de la structure 
lineique a partir des positions sur les deux images d'un ensemble de points de marquage, 
lesdits points de marquage etant ponderes en fonction de leur eloignement par rapport a la 
structure lineique. 

5 

6. Procede de traitement damages selon Pune des revendications 1 a 5, incluant une 
etape de determination d'un ou plusieurs mouvements globaux de la structure lineique, 
lesdits mouvements globaux etant extraits de I'expression mathematique de la deformation 
et une etape de soustraction de ces mouvements globaux a I'expression mathematique de la 

10 deformation appliquee a la structure lineique. 

7. Appareillage de traitement damages, ayant des moyens pour recevoir ou generer des 
images, lesdites images appartenant a une sequence d'au moins deux images presentant 
une surface representative d'un organe ou d'une partie d'organe deformable dans le temps 

15 et appelee surface d'organe, ladite surface incluant des points caracteristiques, notes points 
de marquage, qui se correspondent d'une image a une autre de la sequence, ledit 
appareillage comprenant des moyens de : 

definition d'une structure lineique sur une des images de la sequence, 
calcul des positions des points de marquage sur au moins deux images, successrves 
20 ou non, 

determination des parametres d'une expression mathematique explicite de la 
deformation de Porgane ou de la partie de Porgane observee entre les deux images, a partir 
des positions d'un ensemble de points de marquage sur les deux images et des positions des 
points de la structure lineique, 
25 application de ladite expression mathematique explicite a la structure lineique pour 

definir la forme de la structure lineique apres deformation de Porgane entre les deux images. 

8. Appareillage de traitement damages selon la revendication 7, caracterise en ce que 
ledit organe est marque par modulation spatiale de resonance rnagnetique, ledit marquage 

30 etant visible sur les images sous la forme de lignes de marquage, lesdites lignes de 

marquage se deformant en suivant la deformation de Porgane et etant telles qull existe des 
points d'intersection entre lesdites lignes de marquage, lesdits points dlntersection etant les 
points de marquage.. 

35 9. Appareillage de traitement damages selon Pune des revendications 7 et 8, pour mettre 
en oeuvre un procede selon Tune des revendications 3 a 6. 
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10. Appareillage de traitement d'images selon I'une de revendications 7 a 9 comprenant 
des moyens pour iterer le precede decrit pour deux images successives ou non dans la 
revendication 1, sur I'ensemble des images successives de la sequence d'images. 

11. Appareillage de traitement d'images selon I'une des revendications 7 a 10, 
comprenant une etape de determination d'un ou plusieurs mouvements globaux de la 
structure lineique, lesdits mouvements globaux etant extraits de ('expression mathematique 
de la deformation et une etape de soustraction de ces mouvements globaux a I'expression 
mathematique de la deformation appliquee a la structure lineique. 

12. Appareillage de traitement d'images selon I'une des revendications 7 a 11 comprenant 
des moyens de visualisation de revolution des parametres de la deformation subie par la 
structure lineique au cours de la sequence. 

13. Appareillage de traitement d'images selon I'une des revendications 7 a 12, caracterise 
en ce que ladite structure lineique est une structure lineique definie par I'utilisateur sur une 
des images de la sequence et en ce que la deformation est suivie sur I'ensemble des images 
de la sequence. 

14. Appareil de saisie dlmages de resonance magnetique comprenant : 

- des moyens d'acquisition d'images de resonance magnetique permettant d'obtenir une 
sequence d'images, 

- des moyens de presentation visuelle de ces images, et 

- un appareillage de traitement dlmages selon I'une des revendications 7 a 13. 

15. Produit programme d'ordinateur comprenant des portions/moyens/instructions de 
code de programme pour I'execution des etapes du precede selon I'une des revendications 1 
a 6 lorsque ledit programme est execute sur un ordinateur. 
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10. Appareillage de traitement dlmages selon Tune de revendications 7 a 9 comprenant 
des moyens pour iterer le procede decrit pour deux images successives ou non dans la 
revendication 1, sur I'ensemble des images successives de la sequence dlmages. 

11. Appareillage de traitement dlmages selon I'une des revendications 7 a 10, 
comprenant une etape de determination d'un ou plusieurs mouvements globaux de la 
structure lineique, lesdits mouvements globaux etant extraits de I'expression mathematique 
de la deformation et une etape de soustraction de ces mouvements globaux a I'expression 
mathematique de la deformation appliquee a la structure lineique. 

12. Appareillage de traitement dlmages selon Tune des revendications 7 a 11 comprenant 
des moyens de visualisation de revolution des parametres de la deformation subie par la 
structure lineique au cours de la sequence. 

13. Appareillage de traitement dlmages selon I'une des revendications 7 a 12, caracterise 
en ce que ladite structure lineique est une structure lineique definie par I'utilisateur sur une 
des images de la sequence et en ce que la deformation est suivie sur I'ensemble des images 
de la sequence. 

14. Appareil de saisie d'images de resonance magnetique comprenant : 

des moyens d'acquisition dlmages de resonance magnetique permettant d'obtenir^une 
sequence dlmages, 

des moyens de presentation visuelle de ces images, et 

un appareillage de traitement dlmages selon I'une des revendications 7 a 13. 



1er depot 



1/4 



DEF v/E 

IM(tl) 



LS(tl) 



AP 



MP(tl) 



MP'(tl) 




f(tl/t2) 



LS(t2) 



] \S S IM(t2) 



I 



V^CALC 



v ' MP(t2) 



MP'(t2) 



V^DET 



FIG.l 




FIG. 2 



1 er depot 

2/4 




1 er depot 



3/4 




FIG. 5 



DEV 




FIG. 6a 



LS(t2) 



FIG. 6b 




r a 



8. 



FIG.6c 



FIG.6d 



1 er depot 



4/4 




FIG. 7 



APP 



BIKSTttVT 
NATIONAL BE 
LAPROPftlETI 
INOUSTRICLLC 



BREVET D'INVENTION 
CERTIFICAT D'UTILITE 

Code de la prapriete intetlectijelle - Uvre VI 



N* 11235*02 



DtPARTEMENT DES BREVETS 

26 bis. rue de Saint Petersbourg 
75800 Paris Cedex08 



DESIGNATION DMNVEMTEUR(S) Page N° X. . / J. . 

{Si le demandeur n'est pas I'inventeur ou Tunique inventeur) 



Vos references pour ce dossier 

(facultatif) 


PHFR010003 


N° D'ENREGISTREM ENT NATIONAL 





TfTRE DE L'INVENTION (200 caracteres ou espaces maximum) 

Suivi de la deformation d'une structure lineique definie sur une image d'une sequence d'images d'un organe deformable dans le 
temps. 



LE(S) DEMANDEUR(S) : 

KONINKLI IKE PHILIPS ELECTRONICS N.V. 



DESIGNE(NT) EN TANT QU'INVENTEUR(S) : (Indiquez en haut a droite «Page N° S'il y a plus de trois inventeurs, 

utilisez un formulaire identique et numerotez chaque page en indiquant le nombre total de pages). 



Norn 


ALLOUCHE 


Prenoms 


Cyril f 


Adresse 


Rue 


1 56, Bd Haussmann 


Code postal et ville 


75008 | PARIS 


Societe d'appartenance (facultatif) 


Societe Civile S.P.I.D. 


Nom 




Prenoms 




Adresse 


Rue 




Code postal et ville 


I 


Societe d'appartenance (facultatif) 




Nom 




Prenoms 




Adresse 


Rue 




Code postal et ville 


1 


Societe d'appartenance (facultatif) 




DATE ET SIGNATURE(S) 
DU (DES) DEMANDEUR(S) 
OU DU MANDATAIRE 
(Nom et qualitc du sign ata ire) 

23 Janvier 2001 
C.THOMAS 

Mandataire SPID 422-5/S008 







La loi n D 78-17 du 6 janvier 1978 relative a 1'informatique, aux fichiers et aux libertes s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concemant aupres de I'INPI. 




BLANK PAG£ 



